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Abstract 

A thin anisotropic elastic plate clamped along its lateral side and also supported at a 
small area of one base is considered; the diameter of 6h is of the same order as the plate 
relative thickness ft. <C 1. In addition to the standard Kirchhoff model with the Sobolev point 
condition, a three-dimensional boundary layer is investigated in the vicinity of the support 

which with the help of the derived weighted inequality of Korn’s type, will provide an 
error estimate with the bound cft^/^|lnft|. Ignoring this boundary layer effect reduces the 
precision order down to | lnft|“^/^. 

Keywords: Kirchhoff plate, small support zones, asymptotic analysis, boundary layers, 
weighted Korn inequality. 
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1 Introduction 

1.1 A plate supported over a small area 

Let w be a domain in the plane with a smooth boundary doj and a compact closure tJ = ujVJdui. 
We introduce the cylindrical plate 


U;, = a; X (-ft/2, ft/2) (1.1) 

of a small thickness ft > 0. By rescaling, the characteristic size of oj is reduced to one so that the 
Cartesian coordinates x = (y, z) G x M and all geometric parameters become dimensionless. 
The bases of the plate and its lateral side are given by 

'^h = {x'-y = (ftu y 2 ) Guj, z = ±h/2} , Vh = {x : y G duj, p| < ft/2} (1.2) 
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Figure 1: A plate clamped over the lateral side and a small support area 


respectively. We fix a point O inside uj and place the y-coordinate origin at O. Denoting by 
0 C an open, not necessarily connected, set with a compact closure 6, we assume that the 


plate (1.1) is clamped over the lateral side as well as at the small supporting area 


6h = {x : r] := h G 0, z = -h/2] 


(1.3) 


on the base . Characteristic sizes of 0 are supposed to be of order one, too. In other words, 


the diameter of the supporting zone (1.3) is comparable with the plate thickness h. The bound 
/iQ > 0 for the small parameter h is chosen such that 6^ C for all h G (0,/iq], however if 
necessary, we will reduce /iq but always keep the notation. Throughout this paper we do not 
distinguish in the notation for 9 and 6^ between two-dimensional sets and their immersions in 
along the planes {x \ z = —1/2} and {x : z = —hj"!), respectively. 

The plate is made out of a homogeneous anisotropic elastic material and its deformation 
is caused by volume forces. It is ideally fixed over the set 


T/i = Uft U 9h 


(1.4) 


while the rest of the plate surface, in particular S* = \9hi is traction-free. The clamped 
area is shaded in Figure 

The main goal of this and consequent paper [5] is to examine the influence of the supporting 


area (1.3) on the stress-strain state of the whole plate. To this end, we construct asymptotics 
of elastic fields as h —)• -1-0, prove error estimates, and create a two-dimensional model which 
reflects adequately all principal effects of the small support. We emphasize that clamping along 
the lengthy set Vh C d^lh plays a considerable role in technicalities of our study and the case of a 
traction-free lateral side will be treated in a subsequent paper where, in contrast to the present 
case, the number of small supporting areas and their location become of a major importance. 


1.2 Formulation of the elasticity problem; the Mandel-Voigt notation 

Since the Cartesian coordinate system x = (xi,X 2 ,X 3 ) attached to the plate is fixed, we can 
regard the displacement field u as the column (tti, tt 2 , "Ws) 
u onto the Xj-axis and T stands for transposition. The strain column 


T ■ 

m 


where uj is the projection of 


£ — £ 22 , , £33) 


(1.5) 
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substitutes for the strain tensor of rank 2 with the Cartesian components 


1 f dun 


dui 


and can be computed by the formula 

e{u) = D (V) u 

where V = grad and 


j,k = 1,2,3, 


( 1 . 6 ) 


(1.7) 


di 0 2-1/292 2-1/293 


0 


D {wy = 0 52 2 - 1/291 0 2 - 1/293 


0 

0 


1 dxj 


7^- (1-8) 


0 0 0 2 - 1/291 2 - 1/292 83 

Notice that, according to the Mandel-Voigt notation, the factors 2 I /2 and 2-1/2 introduced 
in (1.5) and (1.8) for the purpose of equalizing the intrinsic norms of the tensor and the column 
of height 6 (cf. [3l|37] and others). 

The stress column 


<7 = (^CTll, (T22, 2>/V2.2‘'V,3,2>/V.3 ) '733) 


T 


of the same structure as in (1.5) is to be found through the Hooke’s law 

a (u) = Ae (u) = AD{V)u 


(1.9) 


where A is the stiffness matrix of size 6x6, symmetric and positive definite. This matrix 
contains elastic moduli, for instance, to an isotropic elastic material there corresponds 



/ 

\ + 2 fj. 

A 

0 

0 

0 

A \ 



A 

A + 2fj, 

0 

0 

0 

A 

A = 


0 

0 

2 /j 

0 

0 

0 


0 

0 

0 

2 fi 

0 

0 



0 

0 

0 

0 

2 /j 

0 


V 

A 

A 

0 

0 

0 

A -|- 2 fi / 

where A > 0 and /r > 0 are the Lame constants. 






In what follows we use matrix, rather than tensor, notation in elasticity known as the Mandel- 
Voigt notation, cf. [3]. In this way we write the equilibrium equations as follows: 

L{V)u{h,x):=D{-V)'^ AD{V)u{h,x) = f{h,x), x £ Qh, (1-11) 


where / = (/i, / 2 , f 3 )~^ is the vector (column) of the volnme (mass) forces. The three-dimensional 
elasticity system (|1.11 ) is supplied with the traction-free boundary condition 


N+{V)u{h,x) := Dies)'^ AD{V)u{h,x) = 0, x E S+, (1.12) 

N~ (V) u {h, x) := D (— 63 )"'' AD (V) u {h, x) = 0, x E S*, 
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where 63 = (0, 0,1)"*" is the unit vector of the outward normal on the bases At the clamped 
parts of the surfaces we write 


u {h, x) = 0, X G 9h, 
u {h, x) = 0, X G Vh- 


(1.13) 

(1.14) 


We further refer to (1.12) and (1.13), (1.14) as the Neumann and Dirichlet conditions respec¬ 
tively. 

The variational statement of problem (1.11)-(1.14) reads: to find a vector function u G 
Hq (Q/i; T/i)^ such that 

{AD (V) u, D (V) x)^^ = (/, , Vx G {Qh-, Tnf , (1.15) 

where ( , is the natural scalar product in the Lebesgue space (Q/i), is a 

subspace of functions in the Sobolev class {Qh) which vanish at set (1.4), and the last 
superscript 3 in (1.15) means that test (vector) function x has three components. 

In view of the Dirichlet conditions (|1.13) and (|1.14[) the Korn inequality |19] 


\v,H^{nh)\\<Kh\\D{v)v-L^{nh)\\, yvGH^{nh-,rh) 


(1.16) 


is valid and the unique solvability of the variational problem (1.15) with any / G (D/j)^ 
follows from the Riesz representation theorem. Although there exist various approaches, see 
da da [mils] and many others, to verify inequality ( |1.16[ ), we still need to clarify the dependence 
of the Korn constant on the small parameter h. To this end, we exhibit in Section several 
variants of weighted anisotropic Korn’s inequalities. 


1.3 Asymptotics and boundary layers 

The Kirchhoff theory of thin plates created more than 150 years ago by means of intuitive 
asymptotic analysis, has got a justification in miscellaneous formulations by various methods 
and at different level of rigor, we refer only to the mathematical monographs p da ED EZl SH] 
and intensive lists of literature therein, although setting aside a vast volume of publications with 
important theoretical and practical results. If an external loading is scaled to provide the elastic 
energy of a plate to gain order 1 = , the energy norm of asymptotic remainders in Kirchhoff’s 

asymptotic formulas becomes 0 {y/h) and this is the best error estimate available within a two- 
dimensional theory of plates. This limitation of the asymptotic accuracy is prescribed exclusively 
by boundary layer effect near the plate edge. Namely, in the vicinity of the lateral side the 
standard plane-stress state in the (yi, y2)-directions, dominant in the midst of the plane, couples 
with a plane-strain state in the (n, 2;)-directions, where n is the normal vector to doj. The latter 
is represented by special elastic fields which slowly vary along the edge, produce strains and 
stresses of order 1 at the lateral side Vh but quickly, at the exponential rate, decay at a distance 
from Vh- We again refer only to the mathematical papers [33 EH] and HDH, where the edge 
boundary layer phenomenon was investigated in elasticity, and to Chapters 15 and 16 of the 
monograph [23] where simplest scalar and general elliptic problems were examined and basic 
principles to construct boundary layers are laid down. 

The energy norm 

{AD{V)u,D{V)v)l{^ (1.17) 
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of the boundary layer equals 0 {y/h) that just predetermines the error bound in an integral 
norm. However, the proximity of the plane-stress state in a weighted Holder norm is much less, 
cf. [311139] . For example, bounds in the pointwise estimates of remainders in two-dimensional 
asymptotic forms of stresses get the same order in h as the main terms and this is known in 
mechanics as “edge effect” in a thin plate. Thus, the error 0{y/h) of the Kirchhoff model is 
mostly due to the intrinsic localization of the edge effect in a c/i-neighborhood of the lateral side 


Vh- 

The two-dimensional structure of boundary layers is kept for a smooth contour du) only. 
Plates with angulate edges were considered in [26] where three-dimensi onal boundary layers 
with a power-law decay rate were detected. However, their energy norm (1.17) becomes o(\//i) 
so that they play a secondary role. 

As it can be easily predicted due to the Sobolev embedding theorem, the Dirichlet condition 
(1.13) leads to the Sobolev condition (3.6) in the limit, that is the overage deflection G 
of the plate vanishes at the point O. However, the small support area (1.3) also provokes a fully 
three-dimensional boundary layer in the vicinity of 9h and provides a perturbation of order 
I ln/i|“^/^ in the two-dimensional model. In other words, the convergence rate 0(| ln/i|“^/^) of 
the rescaled three-dimensional displacements, cf. Theorem]^ to a solution of the limit Dirichlet- 
Sobolev problem (3.7), (3.8), (3.5), (3.6) is very slow and, therefore, unsatisfactory. At the same 
time, the asymptotic structures derived in this paper provide the same accuracy as in the 
Kirchhoff theory. 

In 1^, the same effect of a crucial reduction of the asymptotic accuracy due to a small 
Dirichlet area has been observed and discussed in detail for a scalar homogenization problem in 
a perforated planar domain. 


1.4 The elasticity capacity and self-adjoint extensions 

To the best knowledge of the authors, this paper is the first mathematically rigorous and com¬ 
plete study of the boundary layer effect near a small support area and its influence on the whole 
stress-strain state of a plate. The necessity to examine rapid variations of elastic fields in the 
vicinity of the support area 8 ^ accounts for the following result obtained in the consequential 
paper |5|: the convergence rate to the two-dimensional Kirchhoff solution of the true three- 
dimensional one is extremely slow, of order | ln/i|“^/^, but involving the boundary layer reduces 
the error estimate bound down to lnh\ and makes it acceptable for applications. At the 

same time, the whole boundary layer solution is not known explicitly and can be determined by 
solving a complicated elasticity problem (3.23)-(3.25) in the spatial infinite ply x (-1/2,1/2), 


for example, numerically. In this way a simplification of asymptotic expansions and further mod¬ 
eling take on special significance. 

As usual in the theory of elliptic problems in singularly perturbed domains (see the mono¬ 
graphs and others) the alternation of boundary conditions on a small set leads to 

asymptotic forms engaging singular solutions of limit problems and in our case the Dirichlet 


condition (1.13) on 6 ^ affects the limit problem in uj in two ways. First, it requires to impose 


the point condition w-i{0) = 0, (3.6), for the average deflection w^{y) of the plate. Second, the 


far-held, eligible at a distance from the support 6 ^, contains the Green matrix, cf. Section 3.4 
and see [5j, whose entries are solutions in u with logarithmic singularities at the point y = 0. 
Due to the Sobolev embedding theorem C C in the point condition is well-posed within 
the Hilbert theory, cf. [221 Ch.2] but the Green matrix does not belong to the intrinsic energy 
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space for the Kirchhoff plate a; and therefore breaks the usual variational formulation of the 
model. 

To develop a two-dimensional model of the locally supported plate we will turn in [3] 
to the technique of self-adjoint extensions of differential operators (see the pioneering paper [2] 
and, e.g., the review m)- Namely, we get rid of the complicated three-dimensional structure of 
the boundary layer and simulate its influence on the stress-strain state of the plate by a proper 
choice of singular solutions. In other words, we select a special self-adjoint extension in (w)^ 
of the matrix C (V) of differential operators in the conventional Kirchhoff model described in 


Section 3.2 formula (3.14). Parameters of this extension depend only on the quantity |ln/i| and 


the algebraic characteristics C'^ {A, 0) of the small support area, namely the elastic logarithmic 
capacity matrix introduced heuristically in [25], a 4 x 4-matrix composed by coefficients in 


asymptotics of solutions to the elasticity problem in A, see Section 3.5, and is quite similar to 
the logarithmic capacity in harmonic analysis, cf. |20lllZ| and Remark [l^ and polarization and 
capacity matrices in elasticity, cf. [T] and (Mum. There exist numerical schemes 
evaluate such characteristics. 


to 


1.5 Architecture of the paper 


In Section we derive anisotropic weighted Korn inequality in the plate (1.1) clamped along 
the lateral surface and the small set ( |1.3| ) (Theorem as well as demonstrate that weights 
introduced in the Sobolev norms are optimal. Then we modify our approach to prove Theorem 
[^ which serves for a plate with a traction-free lateral surface but several small support areas. 

The last section is devoted to three-dimensional boundary layers emerging in the vicinity of 
small support areas. First of all, we outline the two-dimensional Kirchhoff model of the plate 


Q.h and the Sobolev point condition 103 ( 0 ) = 0, (3.6) for the average deflection W 3 {y), which 
imitates the entire Dirichlet condition (1.13) on Oh- Moreover, we formulate Theorem from 
[5| on the convergence of the spatial true solution to the Kirchhoff-Sobolev solution in order to 
indicate the extremely slow convergence rate 0(| Inand to emphasize the constitutive 
influence of boundary layers described by an elasticity problem in the layer 


A = m2 X (-1/2,1/2) 


(1.18) 


clamped along the set 6 on its planar base x { — 1/2}. 

To prove the existence of a unique solution in the unbounded domain A but with a finite 
energy, we need a new anisotropic weighted Korn inequality which, however, can be derived by 
means of the same tricks as in Section At the same time, a rigorous derivation of asymptotic 
expansions of the solution at infinity requires both, a formal procedure of dimension reduction 
scheduled in Section 3.2 and a specific application of the Kondratiev theory [17] . cf. 


This 


evident similarity of several approaches in our paper explains its architecture. 


In Section 3.4 we strictly determine the elastic logarithmic capacity matrix which was men¬ 


tioned in Section 1.4 and will be a particular object in [5|. In Theorem 13 we establish its general 
properties. 
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2 Weighted anisotropic inequalities of Korn’s type 

2.1 Some Hardy-type inequalities 

Here below we collect some Hardy-type inequalities that will be used subsequently. 


Proposition 1 Let a{x) he a nonnegative function in (0,1); then the following Hardy inequality 
holds: 


a(x)|u(rE)|^ dx <A 


"l 1 . rl 


a{x) 


a{t)dtj \u'{x)\‘^dx yu G (0,1) with u{0) = 0. 


Proof. Since u(0) = 0 we have 


so that 


u‘^(x) = 2 f u{t)u{t)dt, 

Jo 


j a{x)\u{x)\^ dx < 2 J a{x)(^J \u{t)\\u {t)\dt^ dx. 
Interchanging the order of integration gives for the right-hand side 

2 J |u(t)||u^(t)| ^ J a{x)dx^dt, 


which, by Holder inequality, is majorized by 


a{t)\u{t)\‘^ dt 


iV2r /•! 1 


^JO 


^JO 


a{t) 


a{x) dx\ \u {t)f dt 


il/2 


and this concludes the proof. 


We list below some particular cases of functions a{x), with the corresponding Hardy-type 
inequalities. By taking 

a{x) = 

we obtain for every u G H^{0,1 ) with u(0) = 0 


'0 

which implies the classical Hardy inequality 

rT 


f x ‘^\u{x)\‘^ dx < 4: f {1 — x)‘^\u'{x)f dx, 

Jo Jo 


f x ^|u(x)p(lx<4 f \u\x)fdx 

Jo Jo 

for every u £ H^{0,T) with u(0) = 0. Taking 

a(x) = (1 — ln(l — x)|“^ 
gives for every u £ H^{0,1 ) with u(0) = 0 

f (1 — x)~^\ln{l — x)\~‘^\u{x)f dx < 4 f {1 — x)\u'{x)\‘^ dx, 

Jo Jo 


( 2 . 1 ) 
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which implies the inequality 


pR pit 

/ x~^\h\{x/R)\~‘^\u{x)\‘^ dx < 4: / x\v! {x)\‘^ dx 

Jo Jo 

for every u G H^{0,R) with u{R) = 0. Taking 

a(x) = ln(x/2)|“^ 

gives for every u G 1) with u{0) = 0 


fit 


fi fi 

„-3| ^ /I / „-l| 


/ X '’|ln(x/2)| ^\u{x)\'^ dx < 4 / x ^|ln(x/2)| ^\u {x)\^ dx 
/o Jo 


fit/2 

X ^|ln(x/-R)| ‘^\u{x)\‘^ dx < 4 / x ^|ln(x/i?)| ‘^\u'{x)\^ dx (2.3) 

Jo 


which implies the inequality 

fit/2 

Jo 

for every u G H^{0, R/2) with n(0) = 0. Taking 

a{x) = {x + h)~‘^ with h > 0 
gives for every u G 1) with u{0) = 0 


( 2 . 2 ) 


/ {x + h) ‘^\u{x)fdx<- / {x + h) "^\u {x)fdx, 
lo 9 Jo 


which implies the inequality 

rT 


/ {x + h) ^\u{x)\‘^ dx < - / {x + h) ^|■u (x)f(ix 
/o 9 Jo 


(2.4) 


for every u G H^{0,T) with u{0) = 0. 


2.2 Anisotropic weighted Korn inequality 

In the paper |39] devoted to justification of the Kirchhoff theory of plates, see also the pioneering 
papers [Mj and [8] together with the monographs p uni ED EH HE] etc., it was proved that a 
constant K (w) in the inequality 

|||n;0;,||| < A(w)||Z?(V)n;L2(f2;,)||, yuGH^{nh,Vh) (2.5) 

does not depend on h G (0,1] and, of course, on u, where |||M;n/i||| is the anisotropic Sobolev 
norm 




( \'^yUi\‘^ + h 


. i=l 


dui 


dz 


+ 


dm 


dyi 


+ \Ui\ 


dx. 


( 2 . 6 ) 










Notice that the distribution of coefficients W in (2.6) is optimal, namely one cannot replace h? 
by /i", with a < 2, without losing the independence property of K{uj). 

In the paper m a convenient modification of norm ( |2.6[ ) was suggested, namely the weighted 
anisotropic norm 


\u; n 


/iINo 


I Wy'^i\ + 

L i = l \ 


e 

|2 


dui 

dz 


+ 


dm 

dyi 


1 , 
H— 

4 


(2.7) 




with the weighting function 


+ + ^\m 


Sh {y) = h + dist (y, duj) 


dx 


( 2 . 8 ) 


Notice that in the middle of the plate the weights in (2.6) and (2.7) are equivalent but in the 
vicinity of the lateral side Vh where the Dirichlet condition (1.14) is imposed, the displacements 
Ui, m and all their derivatives get multipliers h~^ and = I respectively. 

Our further justification scheme requires to insert into the norm (|2.7[) some additional weights 


ShM = {h^ + | 2 /|')-^/' (1 + \Hh^ + | 2 /|')|) 


-1 


(2.9) 


that become big in a c/i-neighborhood of the set 6h and, therefore, take the Dirichlet condition 
(1.13) into account. Note that the function Shk is smooth with any h > 0 and in the sequel it is 
proper to put into (2.7) the weight /i + sq {y), tantamount to (2.8). Here, sq stands for a smooth 
in uJ and positive in uj function equivalent to dist(-, So;) in the vicinity of the boundary du. 

Then we will prove the anisotropic weighted Korn inequality in the following theorem. 


Theorem 2 Any displacement field u G {^h) satisfying the Dirichlet conditions (1.13) and 
( 1 . 14 ), meets the weighted anisotropic Korn inequality 

|||u;D,,|||. < K, {uj) \\D{K)u-L^ (D;,)|| , Vu G Hi {Dh^vnAOh) 


where 


\u] Q 


/i-l !!• 


'^h 




i=l 


q2 

2 ^hl 


dUi 

dz 


+ 


+ \dzm\^ + -^S]i2\m\‘^ 

4 


dm 

dyi 

dx. 




2 I 


( 2 . 10 ) 


( 2 . 11 ) 


with a constant K (uj) independent of h £ (0, ho] and u G 77q (D^,; U OhY 


Proof. Due to a completion argument, it suffices to verify inequality (2.10) for any smooth 


function u in which vanishes near Vh and see ( |1.2| ). We use the stretched coordinates 

i = iv, C) = {h~^y, h~^z) (2.12) 
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and consider the circular cylinder Qhr = 
and the radius R is fixed such that 0 C B 


: B^^ X (—/i/2,/i/2) where B^ is the disk {y : \y\ < p} 
^^ 2 - We write 


||V,n; + /i-2 ||n; (Q,.h)||' = h || V^u; (Qfl)||' + h-^ ||u; (Qk)||' 

= /i 11^ ^ u(x); (Qk)II' < c{R,Lo)h ||l?(Vg)n; (Qr)||' (2.13) 

= c{R,oj) \\D{V^)u;L^ (Q/^R)||' < c{R,uj) ||Zl(V)u; L2(00||' • 

Here, we have applied the standard Korn inequality 

(Qr)II' <c{R,oj) \\D{V^)v,L^ (QR)f , Vu G (Qk;0)" 

which holds true due to the Dirichlet condition on 9. 

Let X G (17°° (M) be a reference cut-off function such that 


X (r) = 1 for r < 1/2 and X (^) = 0 for r > 1, 0 < x ^ 1- 

We set Xh{y) = 1 — x {{^R)~^\y\) and, in view of \XyXh (y)| < ch~^, obtain 

\\D (V) (Xhu ); (L!;,)|f < c(||Zl(V)u; L‘^{nh)f 

+ h-^ ||u; {{MhR \ MhR/ 2 ) X {-h/2, h/2)) If) < c \\D{V)u-, L\nh)\f • 


(2.14) 


(2.15) 


Since Shi{y) < ch ^ in QhR, estimates (2.13) and (2.15) mean that we further can treat the 
product XhU only. This product is still denoted by u but we remember its specific property 


u{y,z) = 0 for \y\ < hR/2. 


(2.16) 


Formulas (1.6), (1.14) and integrating by parts yield 


4 / \£i 2 {u)\^ dx = 
J fih 


> 



dui 

dy2 

2 

+ 

du2 

dyi 


dui 

dy2 

2 

+ 

du2 

dyi 


dui 

dy2 

2 

+ 

du2 

dyi 


+ 2 


+ 2 


dui du2 
dy2 dyi ^ 

dui du2 
dy2 dyi 


dx 


dx 


dui 


dyi 


duq 


dy2 


dx. 


Hence, recalling the strains en (u) and £22 (u), we obtain 

\\XyUi-,L^nh)\f < 2||Zl(V)u;L2(L!;,)|f , i = 1,2. 
We have the Friedrichs inequality in w integrated in z x (—/i/2,fi/2), that is 

(i4/i)|f < Coj \\XyUi;L^ (f^fe)|f • 

Let us explain how inequalities (|2.17|)-p.l8 ) provide the estimates 


(L!0|f < c ||L» (V) u; (Ll;,)|f , 


f = 1,2, 


(2.17) 


(2.18) 


(2.19) 
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where 5oi (y) is given by (2.9) at h = 0. To attach the weight s. we rewrite the function Ui in 


the natural curvilinear coordinates n, s while n is the oriented distance to duj, n < 0 in w, and s 


is the arc length along doj. Setting t = h — n and U (t) = Ui (n, s) in (2.1), we observe that the 
Jacobian of the coordinate change y i—)• (n, s) and its inverse are bounded in the J-neighborhood 
Vs of dijj, J > 0 being fixed small. Since |9|„|nj(x)| < \VyUi{x)\, integrating in s G w and 
z G (—h/2,h/2) converts (2.1) into the inequality 


\s^\i-,L^{iujnVs) X {-h/2,h/2))f <c\\VyUi;L‘^{{ujnV5) x {-h/2,h/2))f (2.20) 


because Sh{y) ~ h + |n|. Formula ( |2.18[ ) allows us to replace a; n V5 with a; in (2.20). 

The weight S'oi ought to be inserted in a similar manner. However, we go over to the 
polar coordinate system (r, (p) in the y-plane and multiply Ui with the cut-off function Xu) (y) = 
X (y/Ruj) where is chosen such that C to. It remains to mention that dy = rdrdip and 

||<9,. ixujUi) ; II < Co: Hn*; (B|^) || , 


and to apply inequalities ( |2.2[ ) integrated over {(p,z) G (0,27r) x (—/i/2,/i/2) and (2.18), (2.17). 

To get an accurate information on the derivatives dzUi and XJyU's is a much more complicated 
issue and we apply a trick from m- To this end, we introduce the cut-off func tion Xh{£} = 
X (2z//i) which is null on the bases and observe that, according to definition (2.8), (2.9) of 
the weights, there holds 


^<h^Sh{y) ^SQi{y)<cs, y G w ^ ^ (0>o] • 




( 2 . 21 ) 


Notice that we have reduced our analysis to the case u = 0 in QhR /2 — ^hR /2 ^ ( hj2,hj2). 
Then we proceed as follows: 


4 / k*3(ii)| dx>h^ I XhSh^S'^^ 

J J 



^y^ 

I 2 = XhSh‘^SiiUi-^^dx - 2h? XhdzXhSl^Shm^dx =: + 21^. 


Since \dzXh {z)\ < c^h we have 


<c^\h? 


q-2 q2 

Xh^h '^01 


du. 


dyi 


^ dxj 


\ 1/2 


\ui\^ dx 


1/2 


< + Cri lls^lFolui; (H;,)||' 

where J > 0 is arbitrary and the last norm has been estimated in ( |2.19[ ). Furthermore, 

2 „- 2 o 2 ,2 f 2 ^^3 9 / -2 o 2 \ jh . jh 

^ ^ 9oi)dx-.I^+h, 


4 


' Qy 


both the integrals getting appropriate bounds with simplicity. First, in view of (2.21) 
1 


r 

re; 


< ^4 (ikii 4) ; (^^4|k + 1^33 (rr) ; (^^h)|k) < c \\D (V) n; (H;,)|k • 
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Second, performing differentiation of weight in y we see that, for y G to and h G (0, Hq] 

1 




tot) 


< ch^Sh (y) ^ Soi (y) 


h + so{y) \y\ 


+ A ) ^ (y) ^ 5'oi (y) 


and, thus. 


< c + \\d,U 3 ; (0;,)|p) <C\\D (V) n; (fi;,)|| 


Collecting relations listed above provides the formula 

<c(l + 5-^) \\D{V)u-,L‘^ inh)f-26l’^ 

so that putting 5 = 1/4 leads to the following estimate of the derivatives in question: 
\\s^^ SoidzUt; {^h/2)\f + \\s~h^ Soidu^ / dyi] 

<C\\D{V)u;L^ (0011', 


( 2 . 22 ) 


however only inside the thinner plate ^h /2 = w x {—hl4:, /i/4) where the cut-off function equals 
one. 

We postpone spreading of estimate (2.22) onto the entire plate and conclude with the dis¬ 
placement U3 itself. 

Operating with (2.4) and (2.3) in the same way as with (2.1) and (2.2), we derive from 
estimates (2.22) of Vytta the formula 

< c ||Zl(V)ri; 


(2.23) 


with a weight required in (2.11) but again in a thinner plate. By the way, we become in position 
to lighten weights in (2.22), (2.23) by the replacement S'og t-G- Shq and write 

(2-24) 

Moreover, we now may forget about the artificial property ( |2.16 ) of the displacement field u. 

To improve the obtained estimates, we apply a method of passing anisotropic Korn in¬ 
equalities from one part of a body to another part which was proposed in 1271 and elaborated 
in m- Let Qh be the cube |x : \yi — ?//| < h/2, \z\ < /i/2} with some center (^*^,0). Ex¬ 
tending u by zero from Qh onto the layer {x : \z\ < hj2}, we assume that y^ G uj and set 
= u{y^ + hrj,hQ where ^ = (?7,C) = {v ~ y^) ^ Qi stretched coordi¬ 

nates, cf. (2.12). We introduce the rigid motion matrix of size 3x6 


1 0 0 0 6-6 
d(6=| 0 10-60 6 

0016-6 0 

and make the following decomposition in the unit cube Qi: 

U{i)=U^ {i) + d{i)U^. 


(2.25) 


(2.26) 
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Here, the last term implies a rigid motion generated by the column 


ZY° = d(Q;) ' / d{0'^U{0d^€ 


(2.27) 


where £ Qi : |CI < 1/4} is a half of the cube and d (Q}) is a Gram matrix of size 6x6, 

symmetric and positive definite. 


d(Q'i) = / dicVdiOdC- 


(2.28) 


Owing to definition (2.27), (2.28), the component U'^ meets the orthogonality conditions 

[ d{O^U^{C)d^ = 0€R^ 


(2.29) 


which, as known (see the proof of Theorem 3.3.3 m or Theorem 2.3.3 EZ!), assure the Korn 
inequality on the intact cube 


(Qi' 


< K 


D{V^)U^-,L^ (Qi) 


(2.30) 


Due to the central symmetry of the integration domain in (2.28) the matrix d (Q}) is diagonal. 
Hence, formulas (2.27) and (2.25) immediately show that 

\U^\<cm-L^Q[)\\, i = l,2, \ui\<c\\Us-,L^Q[)\\. (2.31) 

The component is estimated in the following way: 


\uS\=d7^ 


*66 


= d 


-1 

66 


“ dgg^ 


< c 


(0 - (0)d? 




ryf 1\ dZ72 frjl l\dUi , 


9^1 r2,^r 


dr]2 




+ 




drii 


(Q'l) 


Note that integration by parts did not bring a surface integral because [r]f /2 — 1/8)^.=^1/2 ~ 
Referring to the formulas 


^ dC\2 32 

we finally obtain in a similar manner that 

dU 


\uU<c( 


+ 


2 32 


9Uz.r2l^r 


= 0 , 


C=±l/4 


dry, 


(Q/) 


i = 1,2. 
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Now we return to the x-coordinates and the displacement held u. The decomposition (2.26) 


determines the component which, according to (2.30), gets the estimate, cf. (2.13), 

2 




+ h 


-2 


u^;L^{Qh) <c\\D{V)u;L^ 


(2.32) 


Then we calculate 


SifiUi] L 


< c[ h 


,-2 




< cl \\D (V) u; U 


+ Sh {jP) Sih (y^) mes3(Qh) + |^6-i| + |^o | 


^+ \\sJpS3hUi-,L\^j,)f+ h^ 


«-lQ r2/^/ 


(2.33) 


+ h 

< c (^\\D {V) u] L 


.-iQ ^'»3. r2r^> 


+ IIp; 


+ II j T 

/ ll|2' 


/. r2/^/ \||2 


This calculation needs a detailed commentary. We here and further take into account the 
following formulas for weights: 


Sh (y°) Shq (y°) < sup Sh (y) ^ Shq (y) < CSu (y°) Shq (2/°), 
h^'^Shiyy^ Shqiy) < Cq, 9 = 0,1. 


(2.34) 


The hrst inequality in (2.33) was obtained by using decomposition (2.26) and a direct compu¬ 
tation of the norm in T^(Qh)- The factor mes^Qh = was compensated due to the relations 


Uj] L 


2 fin,/ 


= h-3||^,;L2(Q;)|| 


dxk 


= h 


-1 


dW 




dik 


but the last one still leaves the coefficient h? on the square of the L^(<I2)j)-norm of a derivative. 
Then we applied estimates (2.31 )-(2.32). Finally we recalled dehnition ( |2.11 ) of a weighted 
anisotropic norm^ _ 


Augmenting (2.34) with the relation Shi (y) ^ Sh 2 (y) < h ^ in Qh, we continue as follows: 
|2 

2 


\\Sh^S2hU3;y 


< c{ h ^ 

< c 




+ h^Sh (y°) ^ S 2 h (2/°) messQ/j + \ul\^ + 


+ y ||s;,^52hM3;T" 


h? \\sP-SihyyU3]L 


2 an' 


D{v)u^-Lyqh) 

















































where as usual u' = {ui,U 2 )^ ■ Since the left 3 x 3-block of the rigid motion matrix (2.25) is 
annulled by differentiation, treating derivatives of Uj becomes much simpler: 








< C 




+ h^ Kf) 


<c(^||Z)(V)n;L2(Q;,)||' + 




-k 


^^2 j2(^, 

Wi- 




-k \\s,^^ShiVyU3]L‘^{Qh)\\ 


(2.35) 


< c 


Vu^;L2 


2/^/ M|2 


+ h'^Sh {y°) ^ Shi (y°)^ {Pif + 1^51^) 

(p\ ! 2 

\\DiV)u;L\qh)\f + h^\\spShi-^;L\q^\ + h‘^ \\s^^ShiVyU3; 

By the way, one may avoid to present (2.35) down because the (n/i)-norms of dui/dy 2 and 
du 2 /dyi had been estimated in (2.17). However, we observe that the derivatives dui/dyt = £ii{u) 
and dzU 3 = e 33 (u) figure in the stress column (1.5)-(1.7) and collect estimates obtained above 
to arrive at the equality 


|u; 


l<c[\\D{V)u-U 


2 III rtTi' ll |2 

“k 111 '^^) 


Summing these inequalities up over all cubes which are erected from cells of the quadratic net 
of size h in the plane and have nonempty intersection with the plate H/i, yields the estimate 


m: 


nh\\\l<c(^\\D{V)u;L^ (H;,)||' + |||u;H;,/2|||2). 


We combine it with (2.24) getting the result. ■ 

Remark 3 As mentioned above, the optimality of distribution of the weights h and Sh{y) 
in norms ^2.6^ ) and (2.1) is quite known. However, relation (2.13) shows that the norms 
\\uj] (Q/i_r)|| in the small cylinder QhR — X {—h/2,h/2) can be endowed with the big 
factor h~^ but weights in norm \ 2 . 11 ) give the following estimate only: 

\lnh\~^ h~^ \\uj;L‘^ (Q/ii?)|| < c|||u;( 


-q/2 


in 


In other words, it is worth to confirm impossibility of the change Shq{y) i— {h? + \y\ 

( 2 . 11 ) with the simultaneous preservation of estimate ( 2 . 1 (^ . 

Let ijji be smooth nontrivial functions such that ifi (t) = 0 for t ^ (1/2,1). We set Ui (x) = 
V'i (|lnr| / |ln/i|), f = l, 2 , U 3 (x) = 0 and obtain 


/ ^ 1 1 T 

£ii {u;x) = -rr-n^ii 


r |ln h\ \ |ln h\ 


|lnr| 


, i,l = 1 , 2 , £13 (u) = 623 (w) = £33 (w) = 0 , 


where again (t) = 0 for t ^ (1/2,1), that is (|lnr| / |ln/i|) = 0 for r ^ {h,'/h). We thus 
have 




f'/h 


rdr c^nh 


2 \ In h\ 


(2.36) 
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At the same time. 


(/i2 + |y|2)-VV;L2 (n^) 


2 /■^d^ 

/ Inr \ 

= 27rh / 

Jh 

k In /i J 


rdr 
h? + 


> 2'Kh\ In h\ 


'1/2 


*' ([T^^l) I = "''■I '‘I II*' -f-' (5.1)11'. 


{h? + \y\^) + I ln(/i^ + |y|^)|) ^Ui\L?' {Vih) 


rVh 


= 2'Kh 


A 


|lnr| 
|ln h\ 


1 + 


ln(/i2 + \yf) 


-2 rdr ^ c^h 


/j2 _|_ ^2 I/j| 


Glancing over these formulas convinces that logarithms cannot be eliminated in the weighted 


norm (2.11). 


2.3 Traction-free edge of the plate with several support areas 

The approach applied above and described at length in the review paper m helps to derive 
asymptotically exact weighted anisotropic inequalities of Korn’s type without requiring the 
lateral side of the plate to be clamped. Let us outline derivation of such inequalities. 

To determine small clamped area on the lower base we fix some points y^,... ,y'^ inside 
for j / k, and put 


0/j = U • • • U = [x-. rj ■.= \y-y^\ < Rh, z = -h/2} . 


(2.37) 


Real supporting sets, of course, may be bigger, for instance 6j^ D i?-)) like in (1.3), but the 
Dirichlet condition 

u{x) = 0, X ^ Qhi (2.38) 

is sufficient for our purpose. 

Since the resultant inequality is sensitive to the number of supporting sets 9j^, we focus on 
the case 

J > 2. (2.39) 

Notice that the Korn inequality remains the same for J = 2 and the most realistic case J = 3. 
This and the peculiar case J = 1 will be commented in Remark]^ 


Theorem 4 If J > 2 in (2.31), then any displacement field u G Hq satisfies the 

weighted anisotropic Korn inequality 


|||u; QhWle < Kq {to) (1 + |lnh|) ||Z1 {V)u-,L^ (L!;,)|| 
with a constant Kq (lj) independent of h £ (0, ho] and the weighted Sobolev norm 


(2.40) 


r 2 


I'U'j 1 10 — 


'Uh 


i=l 


Y, \VyUif + h^S 


2c!2 

hi 


dUi 


dz 


+ 


du?. 


dyi 


+ sh \uA ) Idzm] 


+ h 


dx, 
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where 


(2.41) 


Shg (y) = max {S’/iq {y -y") : j = 1,..., J} 


and Shq are given in (2.9). 


Proof. We take a smooth displacement field u satisfying (2.38) and by the same means as in 
Section 2.2, see, e.g., (2.13), impose the subsidiary conditions, cf. (2.16), 

u{y,z) = 0 for \y - y^ < hR/2, j = 1,..., J. (2.42) 


First of all, we employ an elegant device from [49] and define the vector function U with com¬ 
ponents 

l^i{y,0 = Ui{y,hC), i = l,2, UsiyX) = hu3{y,hC) (2.43) 


in the vertically inflated plate = {{y, () : y € lo, |C| < 1/2}. The crucial property of (2.43) is 
expressed by the relation 


\D{V)u-,L^ (b?;,)|| 



dui 

^ 1 

dui dus 

2x ^ 

dui du2 

2 

duz 

2 - 

i=l ^ 

dyi 

+ 2 

dz dyi 

) + 2 

dy2 dyi 

+ 

dz 



dydz 


= h 


> h 


/ 

J Hi 


dUi 

2 1 

dUi dUz 

1 

dUi dU2 

2 

dUz ^ 

L i=l ^ 

dyi 

+ -/i 

2 

dC ' dyi 

J + 2 

dy2 5yi 

+ /i ^ 

dz 


dydC (2.44) 


r 2 


L 


E 


2=1 


dUi 


dyi 


1 

+ 2 


dUi ^ dU3 


dC dyi 


1 

+ 2 


dUi dU2 
+ 


dy2 dyi 


+ 


dU3 


dz 


dydC 


= h\\D{Vy,dc)U-,L^ (fii)|| 


The term U'^ in the decomposition, cf. (2.26), 


with the column 


Uiy,C)=ld^ (2/,C) + d(y,C)^f° 

U^ = d (Qi)-^ [ d (y, o’" U (y, C) dydC e 


(2.45) 


meets the orthogonality conditions (2.29) under the replacement Q} i—)• and, therefore, using 

the Hardy and Korn inequalities (2.2) and (2.30) yields 


r-^{l + \lnrj\)-^U^-L^ (Hi) 


-1 7/-L. r2 


< c 

< c 


U^;H^ (Hi) 


D{Vy,d^)U^-,L‘^ (Hi 


(2.46) 


= c\\D{Vy,d^)uj-,L^ (Hl)||^ 


Now in view of (2.45) and (2.42), we write 




d {y, C)"^ {y, C) dydC =: G 


/Q 


(2.47) 


hR/2 
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where ~ {(^jC) • I?/< hR/2, |(^| < 1/2| is a circular cylinder and, according to 

(2.46[), the right-hand side admits the estimate 


\R^\ <c{mes3Qlj^^^y/^h{l + \lnh\) r^. ^ (1 + |lnrj|) ^ (Q^^/a) 

< ch^{l + |lnh|) \\D{Vy,d^)U;L^ (12i)||^ 


(2.48) 


Following [2^, see also [lOl §3.4], we sum up equations (2.47), j = 1,..., J, and obtain the linear 
algebraic system 

M{h)U^ = R := + ■ ■ ■ + e (2.49) 


with the 6 x 6-matrix 


M (h) — d(Qjj^/2) “I + ^{QhR/ 2 ) ■ 


(2.50) 


Each of summands in (2.50) is a Gram matrix, symmetric and positive definite. Moreover, 


by virtue of (2.25) and (2.28), a simple calculation shows that 


d(QL/2) = ( d (y^ 0)' d (y^ 0) + :^t^t + O (h) 


!'hR/2) 


1 


12 


(2.51) 


where O (h) stands for a 6 x 6-matrix of size 6x6 with the natural norm of order h and 

1 


0 0 0 1 0 0 

t= I 0 0 0 0 1 0 

0 0 0 0 0 0 


12 


/- 1/2 


CdC, 


TT 


= meS3 


(2.52) 


< ch 


-2 


(2.53) 


We observe that under requirement (2.39) the inverse of matrix (2.50) satisfies the estimate 

In fact, the 6 X 6-matrix = d (y-^ , O)"*" d (y-^ , O) -|- 3 ^t'''t on the right of (^2.5lL is symmetric 
but only positive. However, in view of ( |2.50[ ), estimate (2.53) is a direct consequence of the 
following fact: 


b G 


and b~^M^b = 0, y = 1,..., J 


= 0 . 


(2.54) 


Moreover, the premise in (2.54) is equivalent to 


t 6 = 0 and d ,0) b = 0, j = 1,..., J. (2.55) 

We put the coordinate origin y = 0 at t he po int y^ and direct the yi-axis through y^ s o that 
yf > 0 and y^ = 0. Owing to (2.52) and (2.25) the equalities t 6 = 0 and d (y^,0) 6 = 0 in (2.55) 
assure that 64 = 65 = 0 and 61 = 62 = ^3 = 0. Furthermore, 

0 0 0 0 0 0 

d{y^,0)-d{y\0) = I 0 0 0 0 0 -yj 

0 0 0 0 y| 0 


and, therefore, be = 0 due to (2.55) with y = 1 and y = 2. So (2.53) is proved. 


From (2.49) together with (2.53) and (2.48) we derive that 

|Z7°| < c(l + |ln/i|) \\D{Vy,d,)K;L‘^ (Oi)|| . 


(2.56) 
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Thus, representation (2.45), estimates (2.46), (2.56) and relations (2.44), (2.43) between and 
u give the following Korn inequality 


< cil + \lnh\y\\DiV)u-,L‘^{nh)\ 


(2.57) 


with the only exception: after returning to x and u the right-hand side of (2.44) gains the 
integral [ { dzU^ l'^ dx i nstead of J\dzU 3 \'^ dx as in (2.6). This gap is filled easily because, 
according to (1.5), (1.8), the expression \\D (V) u; {Qh)\\ on the right of (2.57) includes 
/ |e 33 (u)l^dx = / |(9^U3| dx. 

A repetition of calculations in Section 2.2 (word by word but with s/i = 1, the basic relations 


(2.21), (2.34) and (2.42) being preserved) allows us to introduce the weights Shq given in (2.41) 
into the norm. Note that the real reason to put s/j = 1 is the absence of the Dirichlet condition 


1.14) on the lateral side Vh and the consequent impossibility to apply the Hardy inequalities 


(2.1), (2.4). 


A completion argument again completes the proof. ■ 

Remark 5 The Dirichlet condition (2.38) at small support zones, s ee {2.31 ) and (2.39), eannot 
maintain the Korn inequality without the factor 1 -|- |ln/i| as in (2.4C^ . To corroborate this 
statement, we use the test function Ui{x) = n/=iX (|lii (^j/7?)| / |ln/i|) where x is taken from 
( 2 . 14 ) so that Ui {x) = 1 if all rj = \y — yj] > y/hR, and Ui {x) = 0 if some rj < hR. Clearly, 
for u = CiUi, we obtain 

||u; (H;,)||' = ^ ^ ^2 c > 0 , 

and, similarly to 

" ^ ^ ^ -R^\lnh\^Jh 

The desired inference follows. 


'■2 |ln/i| 


Remark 6 If J = 1 and y^ = 0, the matrix M. (h) = d(Q)j^y 2 ) (2.50) becomes 

TT^h^R^ diag {1,1,1, i (i + h^R^) ,\{\ + h^R^) , \h^R^] 


and gets the right-hand bottom entry O [h'^). In this way (2.53) loses validity and estimate p.56|) 
alters erucially for Uq. Thus, the resultant inequality requires a serious modification, too (cf. ]3^ 
and ® § 3.4 and ^5.2j). We mention the displaeement field u{x) = (l — x {fifu)) 92 , 91 , 
whieh satisfies the relations 

11^. 2,2 > c > 0, /1-3/2 \\D (V) u; (H/,)|| < C 

and indicates a power-law growth of the Korn constant as h ^ -|-0. 
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3 


The convergence theorem and the three-dimensional bound¬ 
ary layer 


3.1 The Kirchhoff model with the Sobolev point condition 

In [5] we will present a detailed procedure of dimension reduction with turns the elasticity 
problem (1.11)-(1.14) into the two-dimensional Kirchhoff model of an anisotropic plate (1.1) 
clamped over the lateral side Vh and the small area 6 h, cf. (1.13) and (1.14). We however 
start with a classical and simple result on convergence which can be achieved by any of known 
methods, cf. monographs [91 uni 1211 EZl SH] and other literature. 

We assume the following representationj^ for the right-hand side in (1.11) : 

fi{h,y,z) = i = 1,2, fs {h,y, z) = h^/'^gsiy) (3.1) 

with g = ( 51 ,5-2, ffs) £ . 

If the rigidity matrix A has the form (1.10) in the isotropic Hooke’s law (1.9), the average 
longitudinal displacements w' = {wi,W 2 )~^ is solution of the two-dimensional elasticity system 


- gAyw'iy) - (A' yS/y w'{y) = g'{y), y £uj. 


with g' = {gi,g 2 )~^ and the deflexion W 3 is solution of the bi-harmonic equation 

(y) = 53 (y), 2 / e w. 


fi X + g 2 


3 X + 2g 


(3.2) 


(3.3) 


Here, Vy = {d/dyi,d/dy 2 y , Ay = is the Laplace operator in the y-variables, and the 


coefficient 


A' = 


2Xg 


(3.4) 


A -|- 2 g 

is computed through the Lame constants A > 0 and g >0. 

The condition (1.13) on the lateral side Vh, see (1.2), requires for the Dirichlet condition on 
the contour du 


Wi{y) = 0, i = 1, 2, W 3 {y) = 0, dnW 3 {y) = 0, y € doj. 


(3.5) 


where dn = and n = {ni,n 2 )~^ is the unit vector of the outward normal at duj. Finally, 


the support area 6 ^, see (1.3) and ( |1.13 ), is reflected by the Sobolev point condition 

W 3 (O) = 0 . 


(3.6) 


It is well known, see, e.g., [6], that for any g G (w)^ the problem ( |3.2| ), ( |3.3[ ), ( |3.5| ), ( |3.6[ ) has a 
unique generalized solution w G Hq (cj)^ x Hq (cu) while wi,W2 G {ui) and W3 G (w \ O) 
but in general W 3 ^ (uj). The next convergence theorem holds and it can be proved by an 

impalpable modification of the standard approach for a plate without the small support area 6 h- 

Theorem 7 The rescaled displacements h^^‘^U 3 {h,y,h() and h}/‘^Ui{h,y,hC,), i = 1,2, in the 
three-dimensional problem . 14 ^ with the right-hand side ( fg. l\) converge in x 

(—1/2,1/2)) as h ^ -|-0 to the functions W 3 {y) and wAy) — C^^r-^iy), i = 1 , 2 , respectively, 

^Ui 

where Q = h~^z is the stretched coordinate and w = {wi,W 2 ,W 3 )^ G Llg (w)^ x Hq {u) is a 
solution of the two-dimensional Dirichlet-Sobolev problem (3.2), {3.3(, (3.5), (3.6). 


^In this assumption will be weakned quite much. 
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The main unfavorable result of [S] reads: the convergence rate in Theorem is unaccept¬ 
ably low, namely 0{\\nh\~^^‘^) and thereafter we prepare for an elaboration of the asymptotic 
structures of elastic fields in O/j near 9^- 


3.2 Sketch of asymptotic expansions in a thin plate. 

Theorem remains valid for an anisotropic plate, i.e., with arbitrary rigidity matrix A in the 
Hooke’s law, where the differential equations (|3.2|) and (3.3) are replaced with 

(3.7) 

(3.8) 


C'{Vy)w'{y) = g'{y), y £ uj, 
^3C^y)w3{y) = gsiy), y £ ui, 


where the 2 x 2 -matrix CJ of second-order operators and the scalar fourth-order operator /I 3 are 
given by 


/:'(v,) = p'(-v,)' . 40 p'(v,), 


C 3 {^y) = {Vy)'^A^V3 (V 


(3.9) 

(3.10) 


Here, the symmetric and positive definite 3 x 3-matrix is computed as follows: 


and 


A ^{yy) ^{yz)^(zz)^{zy)^ ^ 


^'(Vy)=(Q' 02 ’ ^siVy) = (2-^/^dl2-^/^dldid2y . (3.11) 

Calculations of asymptotic expansions and the differential operators (3.9), which are entirely 
adapted to the Mandel-Voigt notation, will be presented in [5] but also can be derived by any 
asymptotic procedure for an asymptotic analysis of thin plates, for example, [35] and [371 §4]. 
In Section 3 we will apply the ordinary asymptotic expansion of the solution to the problem 
(1.11)-(1.14) with the right-hand side (3.1) 


^(ra) ^{yz) 

^(zy) ^(zz) 


u(h, x) ~ ^ ^ ^ 

p=0 


(3.12) 


written in an unusual form, i.e. with the help of the following 3 x 3-matrix differential operators 


0 0 0 


^°(C,Vj,) =000 


1 0 -C5i 
W^((,Vy)=l 0 1 -Cd2 


(3.13) 


0 0 1 / \ 0 0 0 
ir^K.v,) = (-cn3,2V2 (| _ Xj i,jv(v,) 

where I 3 = diag{l, 1,1} and J = diag{2“^/^, 2“^/^, 1} are diagonal matrices and V (Vy) is the 
6 X 3-matrix composed of the blocks (3.11) 


V(Vy) = 


V'(Vy) ©3x1 
©3x2 7?3 (Vy) 
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and the null matrices Opxq of size pxq. Coherently with (3.14), the matrix differential operator 
C{Vy) of the affiliated system (3.7), (3.8) involves the following block-diagonal matrix A of size 
6x6, 

^ = diag{7l0, g^0}, C (V) = diag{£' (V,), £3 (V^,)}. (3.14) 


The operators (3.13) are selected such that 

= 0 , = 0 with 9 = 0 , 1 , 2 . 

It is not possible to annul both and but {(, Vy) is fixed such that 

(Ff, Fiy = a (V,) and F| = 0, = 0. 

Moreover, 

/■ 1/2 

/ FldCFG^^AG^ =L^{\/y). 

4-1/2 


(3.15) 


(3.16) 


Let us hint at the choice of operators (3.13) in (3.12) even if a detailed description of the 
dimension reduction procedure will be given in [5]. The differential operators L (V) and (V) 
on the left in ( |1.11[ ) and (1.12) admit the decompositions 

L(V) = (a^) + h-^L^ (V^, d^) + (Vy), 

/V±(V) = (a^) + {Vy ), 

where C = h~^z G (—1/2,1/2) is the stretched coordinate, d(^ = d/dC, and 

L°(9c) = ^(0,0,-9c)^7lF(0,0,5c), (Vy) = D{-Vy,0)'^AD{Vy,0), 

{Vy, 5c) = D{0, 0, -dc)'^AD{Vy, 0) + D{-Vy, 0^AD{0, 0,5c), 

(5c) = D{±e3)'^ AD{0, 0,5,), {Vy) = D{±e3)'^ AD{Vy, 0). 

Then we have 
3 

L hPWP = h-^L^W^ + h-^ (L°IT^ + L^W^) 

p=0 

+ /i° {lV + + L^W^) + {L°W^ + + L^W^) 

+ {L^W^ + LV) + =: Y 

q=0 

3 

Y + /i° {N^^W^ + iV^=^IT°) 

p=0 

4 

+ h^+ =: Y h'^~^G'^^. 

q=0 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


The differential operators {C,Vy) in (3.12), (3.13) as well as (3.9) in (3.19), (3.20) are defined 
uniquely through the matrices D{V) and A in (1.8) and (1.9). For general elliptic problem, an 
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algebraic procedure to construct asymptotics type (3.12) and the resultant operator >C(Vy) in 
(3.14) is developed in [28], see also m §2 Ch.5] and the review [Mj. This procedure serves for 


2 Ch.5] and the review [3 
constructing asymptotic expansions in thin domains and in unbounded domains with cylindrical 
and layer-shaped outlets to inhnity, cf. El §2] and Section [3^ below. 

Finally it should be mentioned that explicit formulas for {C,S/y) w{y) and higher order 


terms in (3.12) are available but are of no further use. 


3.3 Boundary layer effects 


Although Theoremsingles out the solution w of the limit problem (3.7), (3.8), (3.5), (3.6) in w 
as a limit of the solution u of the original problem in the thin plate fl/i, the Dirichlet conditions 


(1.13) at the small support zone 9h maintain only the point condition (3.6) for the deflection ws 


and are not reflected in the problem for the longitudinal displacements w'. At the same time, the 


components wi and W 2 of w' leave small discrepancies in conditions (1.13) unless accidentally. As 


usual, a variation of boundary conditions on a set with a small diameter brings boundary layer 
effects, see [23( Ch.5] for general apprehension. However, due to the assu med c o mpar ability of 


the plate thickness h and diameters of the boundary layer in problem (1.11)-(1.14) exhibits 


a very specific and intricate structure. To describe it, we introduce the stretched coordinates 


(2.12) and observe that changing x i—^ and putting h = 0 transform the domains and 9^ 
into the layer (1.18) between the planes x {±1/2} and the set 9 C S_ of unit size 


respectively. We also denote S, = S_ \ 9. 

Let us recall the following weighted anisotropic Korn inequality which is proved in m , see 
also |32| and m, and looks quite similar to (|2.40|) and (|2.10|). 


Lemma 8 For any smooth and compactly supported vector function v satisfying the Dirichlet 
condition on 9 the inequality 


\u-,Vo^{A-9)\\<c\\D{V^)v,L^ (A)| 


is valid, where 


\u-,Vo\A-,9)f = 


Li=l V 


dvi 


dC 


+ 


dv?. 


dr]i 


+ \y 



(3.21) 


(3.22) 


+ ± S2 




df, 


_ _ 2 I 

the weighted space Vq (A;0) is a completion of (A \ 9) with respect to the norm (3.22) and 

Skiv) = (1 + (1 ± In (1 ± , P=\r]\- 

Our aim is to investigate the elasticity problem 


AD{^^)v{0 = F{0, 

ee A, 

(3.23) 

(Dies)'^ AD (Vg)u(0 = G+(0, 
\zi(-e3)^AZ)(V5)u(0 = G«(0, 


(3.24) 

^^(0 = 0, fe9, 


(3.25) 
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in layer (1.18) clamped along the area 6] we assume that 


E (A), S^^Gf E (S±) , i = 1, 2, 

/• 1/2 

Fz{v,0 = Fl{^,0 + Fl{r}), / F30(7y,C)dC + E^3(^) = 0> (3.26) 

J-l/2 ± 

E (m2) , F^€L\A), G±EL2(S±), 


where G is an extension of G* over 6 (observe that u = 0 on 0). 

By Lemma we c an give a variational formulation to problem ( |3.23 )-(3.25). 

Multiplying (3.23) scalarly with a test function u E (A \ 0) , smooth and compactly 
supported, and integrating by parts with the help of the Neumann boundary conditions (3.24) 
lead to the integral identity 


{AD (Vg) V, D (Vg) u)^ = {F, u)^ + ^ (G±, n)^± . (3.27) 

± 


Then the following result holds. 


Proposition 9 Under conditions (3.26) , pro blem (3.23)-(3.25) has a unique w eak so lution v E 
Vq [A] 9)^ verifying the integral identity {3.21) and the norm ||u; Vg^(A; 0)||, see (3.22), does not 
exceed the sum of norms of functions in {3.26) multiplied with a constant. 


Proof. By virtue of the orthogonality condition in (3.26), the right-hand side F{u) of (3.27) 
can be rewritten as follows: 


F{u) = ( {Fi,Ui)^ + Y{Gh'>^i)T:± 1 + 


(3.28) 


2 = 1 


± 


+ (F1,U3 - Us) ^ + Y (<^3 > ^3 - W3)s± 


where us {rj) = fY/ 2'^3 iv^OdC and \\S2U3-,L‘^ (M2)|| < ||S' 2 U 3 ;L 2 (A)||. 

Since d(^us = 633 (u) E (A), the one-dimensional Poincare and trace inequalities in the 
interval (— 1 / 2 , 1 / 2 ) 3 ( integrated over the plane M 2 3 rj, provide that 

||u3 - Us] L^(A) II + \\u3-U3-, 7.2 (S=^) II (3.29) 

< c ||5;j {u3 - U 3 ) ;7.2 (A)II = c|| 5 ^U 3 ;L 2 (A)|| < c ||u; Vq (A; 0)11 . 

These together with the weighted trace inequality 

||5iUi; 7.2 (S±) II < c (||5ia^Ui; (A)|| + ||5iu,; (A)||) < c ||u; (A; 0)|| 

demonstrate that, under conditions (3.26), we have the continuous functional ( 3.28[) o n the right 
in (3.27) if test fu nctio ns are taken from the space Vq (A;0)^. Since inequality (3.21) serves the 
left-hand side of (3.27) to be a scalar product in Vq (A;0)^, the Riesz representation theorem 
proves the thesis. ■ 
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All elements of (A) (no condition on 9 is imposed) belong to (A). Any rigid motion, 
except for the rotation ■^162 — ^261 which makes the integral in (3.22) divergent, falls into (A)^ 


because it can be approximated in norm (3.22) by smooth compactly supported vector functions 
(see [30l|32] for details). We introduce the submatrix 


10 0 c 

d#(r?,C)= I 0 1 -C 0 

0 0 r /2 -r?i 


(3.30) 


of rigid motion and notice that the third and sixth columns are excluded from the original matrix 
d{T]X) in (2.25). Exhibiting a general result in [32], the next assertion detect 




(3.31) 


as the main asymptotic term of the solution v G V(/(A;0)^ under certain conditions on the 
right-hand sides. 


Proposition 10 Let the right-hand sides of problem (3.2S)-(3.25) satisfy the following smooth¬ 
ness and decay requirements 


\VPdlF, (0 I + |VPG+ (ry) | + (77) | < Cp^p 

|V(;9^^E30(0| + iv) I + \VPGl (77) I < cp,p 

(0 I < cp,p 


-2+e^ 
-2+e^ 
—3+e 


7 = 1 , 2 , 


for p > Re, 


(3.32) 


where p,q = Q,1,2,, e G (0,1) and Re is such that 9 C . Then a solution v G Vq (A; 9)^ 
of problem (3.23)-(3.25) given in Proposition^ becomes smooth for p > R^ with any Rg > Re 
and verifies the inequalities 

iO - Vi ( 0 ) I < Cpq/9^ (773 if) - (0) I < for p > R% (3.33) 


where is the rigid motion (3.31) with a coefficient column c^ G depending on F, G^, G* 


and other notation is the same as in (3.32). 


Note that the imposed orthogonality condition in (3.26) allowed us in Proposit ion [9 


to reduce 

the decay requirement on F^ and G^ because in the examination of functional (3.28) we used 


inequality (3.29) without referring to the weighted norm (3.22). Furthermore, these components 


can be compensated by a solution of a problem on the interval (—1/2,1/2) 9 C, which inherits 
the decay properties from F^ and G^, G* in (3.32) and thus does not influence main asymptotic 
terms indicated in (3.33). 

The asymptotic form n = -I- n is detected in [32| by means of the dimension reduction. 


cf. Section 3.1 and an application of the Kondratiev theory HZ! (see also im Ch.3 and Ch.6]) 
together with various weighted forms of Korn’s inequality. The paper [32| furnishes complete 
multi-scale decompositions of elastic fields in a layer and we will precise the asymptotic forms 
of Proposition in Section 3.6 
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3.4 The fundamental matrix of the differential operator £ (V^) 

The most easily understood way to compensate for, e.g., the main part of the discrepancy 

hr^l'^Wi (y) ei = (O) e* + for y ^ 6h 


in condition (1.13) is just to solve problem (3.23)-(3.25) with F = 0 and = 0, G* = 0 but with 


the right-hand side —Wi {O) in the Dirichlet condition (3.25) on 6. Since by an appropriate 


extension the inhomogeneity can be passed over to the right-hand sides in equations (3.23), 


( 3.24[ ), Proposition]^ and our further comments on the rigid motion (3.31) demonstrate that the 
unique solution u G (A)^ of the problem is nothing but the constant vector —Wi (O) e* which 
does not decay as p —)• -|-oo and by no means can be accepted as a boundary layer. Thereupon, 
instead of solutions of fered b y Pro positions and we prefer to employ a solution of the 
homogeneous problem (3.23)-(3.25) which, of course, must live outside Vg^ (A;0)^. Asymptotics 
at infinity of elastic fields with power-logarithmic growth in a layer has been prepared in the 
paper |32] and to achieve the goal we need a certain notation only. 

General results in m supply us with the fundamental matrix of size 2 x 2 


4)' (y) = ^' Inr + tp' (p) 


(3.34) 


of the elliptic 2 x 2-matrix (3.9) of the second-order differential operators. Here, (r, (p) is the 
polar coordinate system on the plane 9 y, T' is a numeral nondegenerate symmetric 2 x 2 - 
matrix and ip' is a smooth matrix function on the unit circle S. By its meaning, the fundamental 


matrix (3.34) satisfies the relation 


- [ M'{y,Vy)^’{y)dsy = h G 

J •y 


1)2x2 


where 7 is any smooth closed simple contour enveloping the y-coordinates origin O and 


M'{y,Vy)=V' {n{y)y MV {Vy) 


(3.35) 


is the Neumann condition operator for (3.9) with the unit vector n = (u,i,n 2 ) of the outward 


normal at 7 . Since the matrix M in (3.34) is defined up to a constant summand, it can be fixed 
such that ^ 

[ ^{y)'^Af'{y,Vy)^'{y)dsy=02eM.^^^. (3.36) 


The scalar forth-order elliptic operator (3.10) possesses the fundamental solution $3 (y) 
which we write down in the convenient form 


1 


^3 (y) = r --^slnr Ftps (<p) 


(3.37) 


where Ta / 0 is a number and ^p 3 is a function on S. Notice that (y) = diag {‘h' (y), 4>3 (y)} 


implies the fundamental matrix of the operator L (Vy) in (3.9) 
In the isotropic case (1.10) we have 

f - In r -P /3yf r -2 Pyiy 2 r~‘^ 

^ ^ dvr y (A' -h 2 /i) V I3y2yir~^ - In r -h /3y|r 


1 .T, / \ 3 A -|- 2/i 2 1 

.,-2 ). «>!>(!/) = 
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where /3 = (A' + 3^) (A' + /r), A and /x are the Lame constants and A' is defined in (3.4). 

We now introduce the 3 x 4-matrix 


-P , 4>' (y) 0 0 

ch«(y)= (d«(-V„0)^<f>(yr) = 0 0 

' 0 0 <^liy) 


(3.38) 


where (y) from (3.34) has size 2x2 and 


^3 (y) — (y) = ^ 3 yi in r+ ril;l{ip), i = 1 , 2 . 


(3.39) 


To clarify properties of (3.37) and (3.39), we observe that T>^ (Vy) = 2 {Vy)Vy in (3.11) 

and integrate by parts as follows: 


(£3u;3,U3)r + {M 3 W 3 , ( 1 , Vy)u 3 ) = (u;3,/:3U3)r + (( 1 , Vj^) W 3 , A/' 3 U 3 ) . 


(3.40) 


Here, T is a domain bounded by the contour 7 (actually, we need T = cj or T = By and 
A /3 = (A/i), A/i, A/ 2 ), 


Mo (y, V^) = 2-i/2n (y)^ V' (-V^)^ A3V3 (Vy ), 
M (y, V,) = - 2 -V2e7 V' (n (y))^ A3V3 (Vy ). 


(3.41) 


The fundamental solution (3.37) and its derivatives (3.39) satisfy the relations, with f,/c = 1,2 

- (l,A/o‘h3) =1, - ((1, Vy)yfc,A/3<l)3) = 0, 

' (3 42) 

-(1,W3<1>*)^ = 0, {il,Vy)yk,M3^)^ = 5,,k. 


Remark 11 The simplest way to eheck up (3.42) requires to use the Dirac mass 6 (y) in the 
framework of the theory of distributions. Since £3 (V^) 4'3 (y) = <5 (y) by definition, we apply 


formula (3.40) with either u )3 = 1, U 3 = 4>3 or W 3 = yk, V 3 = 4>| and obtain 

-{ilMy)l,M3^3).y = (l,-5)r = 1, - ((l,V,)yfc,AA3<h^)^ = (yfc,a5/%)r = -^k- 


Other relations in (3.42) as well as (3.35), (3.36) are verified in the same way. 


3.5 The elastic logarithmic capacity 

According to a general procedure in [32] we search for a matrix solution of the homogeneous 


problem (3.23)-(3.25) in the form 


r{C) = [I-X 


2 Re 


p=0 


(3.43) 


where the notation from (3.12) is used and V ^ Vq (A;a;i)^^^ is a remainder to be determined. 
Note that Re is hxed such that 9 C 
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We insert (3.43) into the differential equations (3.23) with F = 0 and take decomposition 


(3.15) into account. If p > Re and 1 — xip/‘^Re) = 1, we obtain for the sum E{r]X) of the 
detached terms in ( |3.43 ) that 

Z)(-+ {L^W^ + 

+ {L^W^ + + L‘^W°) + L‘^W^) 4>** (3.44) 

+ {L^W^ + L^W^) ¥ + 


According to the content of Section 3.2, hrst four items on the right in ( 3.44 ) v anish ; reca ll tha t 
V AV (Vn) (r/) = 0 in th e punctured plane M^\0. Owing to (3.16), (3.13) and (3.34), 

(3.39), the last two terms in (3.44) are of order p~^ and p~'^ respectively. 

In the Neumann boundary conditions ( 3.24[ ) we have 

D (ies)^ AD (V^) “ = 4'“ 


First four items on the right vanish again. Recalling that the eqnation (3.20) appeared as a 
result of calculation (3.17), we detect the following relation which is crucial for our further 
consideration: for p G \ 0 


/■ 1/2 

/-1/2 


(V^, 5^) (C, V^) + L\Vr,)W^ (C, V,) dC 


+ (V,) (± 2 ’ (^) j = 0. 


(3.46) 


Let us compose a problem of type (3.23)-(3.25) for the remainder V in (3.43). Since the cut¬ 
off function 1 — X {‘^p/^o) is null near the ry-coordinates origin where <I>^ becomes singular, the 
right- hand s ides F and G~^, G* in the pro blem are smooth and, outside a big ball, coincide with 
(3.44) and (3.45) respectively. Formula (3.46) furnishes the representation F 3 (^) = F^ + 
F^ (ry) for p > 2Re where 

U/2 


Fi (p) = ej / (V,) (C, V,) (ry) = 0(p-"), 


/-1/2 


F^ (0 = ej (L^ (V„ d^) (C, V,) + (V,) ((, V,)) (p) + 

+ejL^ (V,) (C, V,) 4 .« (ry) - Fi (p) = 0(p-^). 


Note that 


(ry) = (V^) W" (± 1 / 2 , V^) 4>» (ry) = 0(p-^}, p > 2Re 


Thus, all hypotheses in Proposition]^ are fnlfilled and then remainder V ^ (A;a;i' 


3 x 4 


(3.47) 


(3.48) 


exists. 


Formulas (3.47) and (3.48) also allow us to derive from Proposition 10 the representation 

r{C) = dHp,C)c^ + rio (3.49) 
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where the rows Vi, i = 1 , 2 , and V^ of the 3 x 4-matrix V meet the inequalities 


(e) I < Cp,p% (0 I < for p > Rl 


(3.50) 


with any e G (0,1). 

We call the numerical matrix C'^ = {A, 6) of size 4x4 the elastic logarithmic capacity 

(matrix), cf. [25j for the isotropic case. It depends on the stifness matrix A in the Hooke 
law jm and the shape of the clamped zone 6 and is defined uniquely through decompositions 
(3.43), (3.49) of the elastic logarithmic potential matrix V (^) of size 3x4. 

Remark 12 The names introduced above come by analogy with the logarithmic capacity poten¬ 
tial P (p) in harmonic analysis, cf. 12(A \4^ . The function P is harmonic in \ d, vanishes at 
the boundary dd of a compact set d and admits the representation 

P{p) = (inp"^ + C\og{d)) +0{p~^), p^ +00, 

while — In \p\ is the fundamental solution of the Laplaeian in the plane and C\og (d) is 

called the logarithmic capacity. Clearly, P also solves the mixed boundary value problem in the 
perforated layer & = : p £ \ d, |C| < 1 / 2 }, namely 

-AgP(e) = 0, eee, P{i) = Q, f£ddx{-l/2,l/2), ±d^P {p,±1/2) = Q, p£^^\d 


which really looks quite similar to a scalar version of the elasticity problem (3.23}-(3.25} under 
consideration. Note that, as distinct from the standard harmonic capacity in dimension n > 3 
which always is positive, the logarithmic capacity Ciog {d) can be positive (the obstacle d is big) 
or negative (d is small). 

The elastic logarithmic capacity matrix C^, in general, is neither positive, nor negative but 
still symmetric. The latter is proved in j25| for an isotropic layer A with a defect and we here 
serve for a much more complicated anisotropic case. 

Theorem 13 The elastic logarithmic capacity 4 x 4-matrix {A, 9) is symmetric. 

Proof. We insert on both positions in the Green formula for the operator L (V^) in the truncated 
layer A (T) = G A : p < T} and let T —)• + 00 . Since V verifies the homogeneous problem 
(3.23)-(3.25), we are left with an integral over the truncation surface St x (—1/2,1/2), the 
lateral surface of the circular cylinder A (T). We have 


0 = 


/St —1/2 


D{p-\,^)~^ AD{V^)V{i) - {D{n{p),Q)~^AD{y^)V{i)Y dfds. 


(3.51) 

where n {p) = p ^p is the unit vector of the outward normal and dsr^ = pdp is the arc element 
on the circle St- 

We now extract from the integrand in (3.51) all terms of order p~^ which contrib ute to the 
limit. Infinitesimal terms o (p~^), in particular ones generated by the remainder V in (3.49) can 
be removed from the integrand with a clear reason while terms growing as p —)• +00 vanish all 
together after integration in ((p, C) £ (0,27r) x (—1/2,1/2) because the limit does exist. In this 
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way, the integral in is approximated by the expression J (rj) — J where, in view of the 
obvious relation D (Vg) = 0, we have 


/*l /2 / ^ ^ 

Jiv)= (Y,W^C,Vr,)<^Hv) + dHvX)Cn d(-,o] AD{V^)Y,WP{C,V^)^Hv)dC. 

J-i/2 Vg = 0 ^ \P J p=0 


(3.52) 


Similarly to (3.44) and (3.45), in the formula 


D (V^) + {D(;W^ + Dr,W°) + {D^W^ + DrjW^) 

p=0 

+ (D^W^ + n^W^) 


where = D (0, 0, 9^) and D^j = D (V,,, 0), cf. (3.16), the first two terms on the right vanish. 
According to ( 3.25| ) and ( 3.34[ ), (3.39), we have 

{D^W^+P (C, V,) + D^WP (C, V^))$« (7?) = 0{p-P), p=l,2 Dr,W^ (C, V^) 4>« (r/) = 0{p-^). 

Taking ( |3.30[ ) and ( |3.13| ) into account, we write the relations 

and 

D (n ( 7 ?), 0 ) d' (C, V^) = (n (t?) , C, V^) + {n {r,)) 
with the matrices (n) = (dl(C) O)"*^ ^3 of size 6x3 and 


C, V.,) = ( ^ (") • <") , <i‘ (C, V,) = ( 0 1 -fop,] 


Comparing decay rates of remaining multipliers shows that we need to keep in (3.52) the following 
two terms of order p~^ (1 + |lnp|): 


ri/2 

1 - 1/2 


(77(7/), C, v^) (<h«(7/) + S (7/, 0) C^y A {D^W^ (C, V,) + Dr,W^ (C, V^)) dC<^^{p) 


(3.53) 


and 


rl/2 

1 - 1/2 


D 


(77 (7?)) ($« (7?) + d« (7/, 0 ) C«)^ A {D^W^ (C, V^) + (C, V,)) (d) • (3.54) 


By formulas (3.19), (3.20) for C' (Vy), £3 (Vy) and (3.35), (3.41) for J\f', Mj we recall the 
calculations (3.17), (3.18) and conclude that (3.53) is equal to the sum 


'+(d«(7/,0))'c«] AA'(7?,V,) 


+ V 


(#«o))'+4o.0)c«)()0j;;;^j)<fO-)) 


(3.55) 


30 
































where and are submatrices with eliminated third lines. A similar argument converts 
( 3.54[ ) into 

CD# (r?) + 4 (r?, 0) ^ AAo (r/, V,) (t?) . (3.56) 


Now we see that 


0 = i'n) - J ivV) dsr^, 

r-)-+oojg \ / 


(3.57) 


i.e. the limit of the integral on the right-hand side of (3.51), is a linear combination of scalar 


products, listed in (3.36) and (3.42), and their conjugates. Thus, a part of the integral (3.57) 


which involves columns of the matrix 4?^^, vanishes due to ’’orthogonality” conditions in (3.36) 


and (3.42) with 0 on the right-hand side. The other part of the integral which involves columns 
of the matrices S and ‘h**, converts into the difference — CK Indeed, we had reduced 

expression (3.52) to the sum of (3.55) and (3.56) so that after integration in G S it is sufficient 


to apply the ”bi-orthogonality” conditions (3.36) and ( |3.42[ ) with 1 on the right-hand side. 
Since according to (3.51) the integral in (3.57) is null, the theorem is proved. 


3.6 Asymptotics at infinity of elastic fields in layer-shaped domains 


If the right-hand sides F and G~^, G* in problem (3.23)-(3.25) are smooth and decay expo¬ 
nentially as p — 7 - -boo, for example have compact supports, then results in |32j serve for an 
asymptotic expansion of the solution v with a remainder of any given power-law decay O {p~^) ■ 
Hence, we can make the decomposition of the elastic logarithmic potential V a bit more precise, 
namely 


V (0 = (1 - Xi2p/Re)) I E WP (C, V,) ch# (p) + S (p, 0) + tHp)]+V (0 • (3.58) 

\P=0 


The first two terms, of course, are the same as in (3.43) and (3.49), but the new remainder gets 


the faster decay properties, cf. (3.50), 




— ^pqP 


-1+e 




^ ^pq P ■) P ^ 2 Rq . 


The latter is caused by the additional term ( 77 ) in (3.58), a 3 x 4-matrix with the rows 


iv) = P {f) + P ^ InpTy {ip), i = 1, 2, 
T“ ip) = (p) + Inp Tf {p) + (lnp )2 Tf (p ). 


(3.59) 


Let us explain where the lower-order asymptotic terms (3.59) appear from. 


The dimension reduction procedure, quite similar to Section 3.1, leads to the following system 
for the third term in the asymptotic ansatz (3.58): 


P(-V^)^AP(V,,)T#(r?) =.F(7?), pGR^\0. 


(3.60) 


The right-hand side F = {Fi,F 2 ,Fs)~^ is constructed in the same way as (3.17), (3.18) and is 
necessary to compensate for discrepancies of the expression H, see a comment to (3.43), in the 


homogeneous equations (3.23) and (3.24), that are the hfth and sixth terms in (3.44) and the 
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fifth term in (3.45) respectively. In this way, representations (3.13), (3.38) and integration by 


parts in the variable C,, cf. calculation in Section yield 
ri/2 


:F,{r,) = ejDj ' A{D^wHC,Vr,) + Dr^W^ {C,Vr,)) dC^^^irj) = 


1 - 1/2 

. 1/2 


pl/Z 

FM = ejDj AD^W^ (C, V^) 

4-1/2 


(3.61) 


+E 


d 


.1/2 


^ dri. 


D„ 


1 - 1/2 


[D^W^ (C, v^) + (C, V,)) dC^^ir^) = p-^F^i^). 


It is important that logarithms figuring in the matrix (p) due to (3.34), (3.39) and (3.38) 


as coefficients on linear functions in pi, are eliminated in (3.61) by differentiating sufficiently 


many times so that functions (3.61) become positive homogeneous in p of degree —3 and —4 


respectively. However, the logarithm Inp returns into components (3.59) of the solution (p) 


by virtue of the Kondratiev theorem on asymptotics (see |17j and, e.g., |41l Thm. 3.1.4]). To 


confirm this, we observe system (3.60) with F = 0 and find out its vector solutions in the form 
of derivatives of columns in matrix (|3.39[), namely 


E 

2 = 1 


{a\di + 0 * 202 ) ^ iv) 
(a?+*0i+a2+*02) iv) 


(3.62) 


As shown in |4H § 5.4], a ny solution with the same positive homogeneity degrees —1 and 0 as in 
(3.59) takes form (3.62). Furthermore, in the cases o^ 7 ^ 0 and of 7 ^ 0 the third component of 
(3.62) stays linearly depe nden t on In p. Hence, the above- ment ioned theore m pres cribes to search 
for a solution of system (3.60) with the right-hand sides (3.61) in the form (3.59) while tJ**, 
and may vanish only under some orthogonality conditions in (S)^ for the angular part 
F^ (99) = (T")* {if) ,F 2 (99) , 0^3 (</’))"'" of F {p). We do not examine these conditions and specify 
(3.59) because the component T is indicated in (3.58) with an auxiliary technical reason only and 


will be excluded from the final asymptotic formulas for the solution u of problem ( |l.ll[ )"( [l.l4 ) 
in the paper [5] 


We emphasize that a solution of system (3.60) is defined up to a summand of type (3.62). 


A unique in decomposition (3.58) of the elastic logarithmic potential V depends on the whole 


data of the problem, in particular on the clamped area 6, and is specified by means of theorem 
on asymptotics in layer-shaped domains, the most challenging assertion in the paper 
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